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Very recently a non-trivial magnetization plateau at 1/4 of the satura- 
tion magnetization was observed in the S = 1 spin ladder BIP-TENO. In 
our previous work we proposed a possible mechanism of the plateau based 
on the second- and third-neighbor exchange couplings which lead to frus- 
tration. In order to confirm the realization of the mechanism, we compare 
the temperature dependence of the magnetic susceptibility and some criti- 
cal magnetic fields obtained by the numerical calculation for the proposed 
model with the experimental results. 
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1 Introduction 



A recent synthesized organic 5 = 1 spin ladder, 3,3',5,5'-tetrakis(A^-iert- 
butylaminoxyl) biphenyl, abbreviated as BIP-TENO Q, is one of interest- 
ing strongly correlated electron systems. It exhibits a field-induced spin gap 
which is observed as a plateau in the magnetization curve. The high-field 
measurement indicated that the plateau appears at 1 /4 of the saturation 
moment. Such a magnetization plateau is predicted in various systems 
A general condition of the quantization of the magnetization was derived 
from the Lieb-Schultz-Mattis (LSM) theorem for low-dimensional mag- 
nets |]l0|. The necessary condition of the plateau is was described as 

Q{S — m) = integer (1) 

where Q is the spatial period of the ground state measured by the unit cell. 
S and m are the total spin and the magnetization per unit cell, respectively. 
Applying this theorem to the BIP-TENO, the 1/4 plateau is the case of 
S = 2 and m = 1/2. Therefore, a spontaneous breaking of the translational 
symmetry(maybe Q = 2) must occur at the plateau. In the previous work 
[11 1 by the present authors two mechanisms of the 1/4 plateau of the S = 1 



spin ladder were proposed, based on the frustrated interactions. In the 
next section we briefly review the mechanisms and show the phase diagrams 
obtained by the level spectroscopy analysis [|l^. The main purpose of this 
paper is to consider the realization of the mechanism at the 1 /4 plateau of 
BIP-TENO, with some quantitative analyses on the critical magnetic fields 
and the temperature dependence of the susceptibility. 

2 Mechanisms of 1/4 plateau 

As an origin of the 1/4 magnetization plateau in the S = 1 spin ladder, we 
introduce the second and third-neighbor exchange interactions. The model 
is described by the Heisenberg Hamiltonian 

H = Ho + Hz (2) 

L L 

i i 
L 

+ J2'^^{Sl4 ■ 824+1 + 824 ■ S14+1) 
i 

L 

+J'i'^^{Si4 ■ 814+2 + 824 ■ S24+2) (3) 

i 
L 

Hz = -HJ2{SI, + SI,), (4) 
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Figure 1: Spin ladder with second and third exchange interactions. 



where Ji,J±, J2 and J3 denote the coupling constants of the leg, rung and 
second- (diagonal) and third- exchange interactions, respectively (Fig. 1). 

Hereafter we put J±=l. Hz is the Zeeman term where H denotes the 
magnetic field along the z-axis and the eigenvalue M of the conserved quan- 
tity J2i ('S'li + 'S'U is a good quantum number. The macroscopic magne- 
tization is represented by m = M/L. In this definition the 1/4 of the 
saturation magnetization corresponds to m = 1/2. In order to explain the 
mechanism of the plateau at m = 1/2, we use the degenerate perturbation 
theory around the strong rung coupling limit Ji,J2,J3 *C 1 [0, 0] . For 
the two spins at each rung, we take only two dominant states; the singlet 
\I'o,o = (I Ti) — I it) and the triplet ^'i^i = | ]]). We introduce a pseudo spin 
T for each rung coupling and map the two original sates singlet ^'9,0 and 
triplet ^1^1 of the S picture to the | -IJ.) and | f|~) states of T, respectively. 
Effective Hamiltonian in pseudo spin can be written as follows: 



eff 



8(Ji - J2) 



EiTt ■ Tt^,) + ^ j:iTr ■ Tr^2 + Ty ■ rr^,) 



(5) 



This is the Hamiltonian of the T = 1/2 XX Z chain with the second- neighbor 
interaction. The magnetization m = 1/2 of the original system corresponds 
to m = in the pseudo-spin system. Referring the well-known features of 
the 5 = 1/2 frustrated XX Z chain [12|, sufficiently large J2 and J3 lead 
to the Neel order and dimerization of the pseudo spins, respectively. They 
correspond to the field- induced spin gap at m = 1/2, that is the 1/4 plateau 
in the original system. The boundary between the spin-fluid and plateau 
phases is of the Kosterlitz-Thouless (KT) type [15|. Therefore, the pseudo- 



spin picture gives two different mechanisms of the plateau; Neel order and 
dimerization, denoted as plateaux A and B, respectively. 

The KT phase boundary can be determined precisely, using the recently 
developed level spectroscopy applied to the low-lying energy levels of the 
finite chains obtained by the numerical diagonalization. We show only thus- 
obtained phase diagrams in the J2-J1 {J3 = 0) and Ji,-Ji{J2 = 0) planes in 
Figs. 2 and 3, respectively. The two gapless phases in Fig. 2 correspond to 
two different ordered states in the classical limit, depending on whether the 
rung or diagonal interactions are dominant Note that there is an upper 
bound of Ji(~ 0.7) for the plateau A, while no bound for the plateau B. 
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Figure 2: Phase diagram on the J2-J1 plane at m = 1/2 
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Figure 3: Phase diagram on the J3-J1 plane at m = 1/2. 
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3 Comparison with Experiment 



Based on the obtained phase diagrams, we discuss the reahstic mechanism 
of the 1/4 plateau of BIP-TENO. The ratio Ji/J_l of BIP-TENO was es- 
timated as Ji/J_L ~ 1.2, fitting the observed temperature dependence of 
the susceptibility x to the numerical calculation for the 5 = 1 simple spin 
ladder Fig. 2 suggests that there is no chance of the Neel plateau for 
Ji/J_L 1.2. In addition the required value of the J2/J1 for the plateau 
is about 0.69 even for Ji/J±_ < 0.7. Ji/J_l ~ 0.69 is too large for the 
realization. Thus the Neel mechanism due to J2 should be discarded. 

Next we consider the possibility of the dimer plateau due to J3. Accord- 
ing to Fig. 3, the plateau would appear if J3/J1 > 0.39 for Ji/J_l ~ 1.2. 
J3/J1 ~ 0.39 is not so far from the realization, because the lattice spac- 
ing along the leg is much smaller than the rung in the crystal structure of 
BIP-TENO Q. Thus we examine the dimer mechanism due to J3 more 
quantitatively. 
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Figure 4: The ratio of the critical magnetic field {B2/B1). 

Using the numerical diagonalization of finite clusters up to L = 8 (16 
spins) and the size scaling technique on the model (2), we estimate the ratio 
of the critical magnetic field Bi that spin gap disappears and B2 that plateau 
begin to appear. The results for J3/Ji=0.4 and 0.5 are plotted versus Ji 
in Fig. 4, gathered with the experimental result of BIP-TENO. It suggests 
that the calculated B2/B1 is closer to that of BIP-TENO for Ji/J± ~ 
1.7, rather than for Ji/J± ~ 1.2. However the experimental estimation 
Ji/J± ~ 1.2 is not so conclusive, because the fitted curve was not obtained 
by the numerical diagonalization of the 5 = 1 ladder, but by some mean 
field approximation for the rung interaction. Thus it would be important 
to fit the direct numerical calculation for the 5 = 1 ladder including J3 to 
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observed x foi' Ji/ J± = We performed the finite-temperature Lanczos 



method |17] to calculate the temperature dependence of x the system 
L = 8. The results for Ji/ J±_ = 1.7 and various values of J3/J1 are shown in 
Fig. 5, together with the experimental result of BIP-TENO. The numerical 
curve well agrees with the measured one for J^/ Ji={)A or 0.5. Therefore, 
the dimer plateau is expected to realize, if Ji/J_l 1.7 and J3/J1 ~0.4 
or 0.5 are satisfied in BIP-TENO. A little difference in B2/B1 between the 
model calculation and the experiment is possibly due to the finite size effect, 
because we used only the values for L = 4 and 8 to estimate it. Indeed the 
ratio B2/B1 of the finite systems is revealed to increase with increasing L 
toward the experimental result. 
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Figure 5: The temperature dependence of the susceptibility x- 



4 summary 

We proposed two mechanisms of the 1 /4 magnetization plateau in the 5 = 1 
frustrated spin ladder. They are described by the Neel order and dimeriza- 
tion of the pseudo-spin system. Comparing the ratio of the two critical fields 
B2/B1 and the temperature dependence of x between the model calculation 
and the experiment, we conclude that the dimer mechanism due to the third 
neighbor interaction is more suitable for BIP-TENO. 
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